Linewidth of single photon transitions in Mni2-acetate 
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We use time-domain terahertz spectroscopy to measure the position and linewidth of single photon 
transitions in Mni2-acetate. This linewidth is compared to the linewidth measured in tunneling 
experiments. We conclude that local magnetic fields (due to dipole or hyperfine interactions) cannot 
be responsible for the observed linewidth, and suggest that the linewidth is due to variations in 
the anisotropy constants for different clusters. We also calculate a lower limit on the dipole field 
distribution that would be expected due to random orientations of clusters and find that collective 
effects must narrow this distribution in tunneling measurements. 

PACS numbers: 76.30.-v,75.45.+j,78.30.-j 



Mni2-acetate ([Mni20i2(CH3COO)i6(H20)4] 
2 CH3COOH-4H20) is a member of a class of high-spin 
molecular clusters that have been shown to cpchibit 
quantum tunneling of the magnetic moment .Era It 
consists of a core of twelve manganese ions with spins 
tightly coupled via superexchange through twelve oxide 
ions, with a ground-state spin S = 10. These clusters 
are separated by acetate and water groups and arranged 
in a tetragonal body-centered lattice so that the nearest 
neighbor distance is 13.7 A and the shortest distance 
between manganese ions in neighboring clusters is 7 A. 
Since it is fairly unusual for a system this large to display 
quantum mechanical properties, it has been the subject 
of much investigation. 

One particularly interesting area of investigation is the 
interaction of the spins with their environment, as re- 
flected in the linewidth of the energy levels. It is not yet 
clear why in some cases the measured linewidth of transi- 
tions provides information about the intrinsic properties 
of the clusters (homogeneous broadening), while in others 
it seems to be due to variations in the local environments 
of clusters (heterogeneous broadening). A comparison of 
the lincwidths we measure for intrawell transitions with 
the linewidths measured in tunneling between wells can 
help to determine the answer. 

The Hamiltonian for the spin clusters is approximately 
givenhyn=-aS'^,-PSi+-f{Sl+Si)-gfiBS-H, where a 
= 0.38 cm-i |nn|j2xlO-4 cm-\ 7 - ±2x lO''^ cm"!, 
and g ~ 2.ElDBl3Q'EfBEj In zero field, states with equal 
\m\ are degenerate. The ground states m = ±10 are 
separated by a barrier of approximately 66 K. 

In this paper, we describe an experiment that measures 
the linewidth of the intrawell m = 10 — > 9 transition (and 



the -10 — > -9 transition) using time-domain terahertz 
spectroscopy. The measurements were made on a pellet 
pressed from small unaligned crystals of Mito-acetate 
prepared according to the procedure of T. LisJlil 

In terahertz time-domain spectroscopy, a nearly single- 
cycle electromagnetic pulse with a length of a few picosec- 
onds is produced when an optical pulse from a mode- 
locked titanium sapphire laser is incident on a micro- 
fabricated photoconductive antenna. This electromag- 
netic pulse is guided through the sample quasi-optically 
and then detected by a second photoconductive antenna. 
This detector antenna is gated by a time-delayed laser 
pulse split from the same pulse that triggered the gener- 
ator antenna. The time dependence of the transmitted 
electromagnetic pulse is measured by adjusting the delay 
between the laser pulses incident on the generator and de- 
tector. The transmitted electric field Et{t) is then Fourier 
transformed to yield Et{uj), the complex frequency de- 
pendence of the transmitted electromagnetic pulse. This 
transmitted electric field can be normalized by the field 
measured with the sample removed from the beam path, 
yielding the complex transmission coefficient, i(a;). 

Generally, terahertz spectroscopy is used to character- 
ize a material response that varies slowly with frequency, 
for which the frequency resolution is not required to be 
higher than tens of GHz. To our knowledge, it has never 
been used to examine linewidths of transitions. This is 
because higher frequency resolution can only be obtained 
with a longer time delay between the generator and de- 
tector laser pulses. It is difficult to align the translation 
stage well enough to allow constant detection sensitivity 
over a large delay range, since the detector is very sen- 
sitive to the position of the focussed laser spot on the 
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FIG. 1: Transmission magnitude at 3 K divided by the trans- 
mission magnitude at 15 K. The peaks are due to transitions 
that are observed at 15 K but are not observed at 3 K because 
nearly all the clusters are in their ground states. 

photoconductor. Through a careful ahgnment process 
we were able to maintain the detector sensitivity within 
a few percent over the entire delay range. It is also diffi- 
cult to eliminate stray reflections from surfaces of optical 
elements and ends of transmission lines. These stray re- 
flections lead to interference fringes in the spectrum. The 
fringes make it impossible to fit the lineshape to a partic- 
ular form (Gaussian or Lorentzian), but it is still possible 
to estimate the linewidth. 

In Figure 1 we show the magnitude ratio of the trans- 
mission spectra for the Mni2-acetate pellet at 15 K and 
at 3 K. This ratio makes the changes in transmission 
due to temperature more apparent. At 3 K, almost all 
the spins are in the ground state (to — ±10), so only 
the transitions from to — ±10 to ±9 are observed. At 
15 K, some higher energy states are thermally occupied 
so more transitions are observed. These absorptions are 
seen as peaks in the transmission ratio. The location 
of these transitions can be fit to yield the first two pa- 
rameters in the Hamiltonian: a = 0.382(4) cm^^ and 
(3 = 8.0(2) * 10"'^ cm~^. These parameters are in good 
agreement with those cited earlier. 

In Figure 2 we focus on the absorption at 300.6 GHz, 
which corresponds to the transition from m = ±10 to 
±9. (The absorption from to = ±9 to ±8 has a similar 
linewidth.) We plot the index of refraction as a func- 
tion of frequency near this absorption at temperature 
T ~ 2.1 K. The index of refraction was calculated di- 
rectly from the (complex) transmission spectrum using 
the following equation. The transmission through a slab 
of thickness = 1.4 mm with complex index of refraction 
n is:E3 

- An 6'""''/'= 

" (^TTFl~~(^?^Pe2«;^' 

This equation assumes that all reflections from the in- 
terfaces are included, even though THz spectroscopy only 
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FIG. 2: Real and imaginary parts of the index of Mni2- 
acetate. The index is calculated directly from the transmis- 
sion using Eq. (0) . The hue is a fit to Eq. (^ . 

collects data within a finite time window after the main 
pulse. However, in this case the reflectivity is low enough 
and the absorption is strong enough that the limited time 
window does not introduce significant error. Notice that 
since this is a magnetic transition, n must be correctly de- 
fined as i/^e/zuoeo, rather than approximated as i/e/eo. 
We stress that our measurement yields the complex index 
of refraction without any modeling of either the lineshape 
or the high and low frequency extrapolations of the re- 
sponse functions. 

We model the absorption using the standard form for 
a magnetic dipole resonance: 

fi/fio = l + — (2) 

We assume that e is constant over this frequency range. 
The curve in Figure 2 is a fit to this form, with = 
e/eo = 1.72, ujo ^ 300.6 GHz, F = 5.5 GHz, and F = 
0.011. The value of F for oriented crystals would be much 
higher, since the electromagnetic pulse only stimulates 
transitions in crystals that are oriented with their c-axes 
parallel to the direction of propagation. 

Although we have used a standard Lorentzian line- 
shape with no inhomogeneous broadening to fit these 
data, it is clear that an equal or better fit could 
be achieved by considering inhomogeneous broadening. 
While we can confidently measure the linewidth, we can- 
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not draw any conclusions about the lineshape. How- 
ever, there has been anothcK-, measurement of this ab- 
sorption line by Mukhin et alB Their data strongly favor 
a Gaussian fit to the transmissionj-with a full width at 
half maximum (FWHM) of 7 GHz^ (For this material, 
the FWHM of ln{\t\), Im{n), and Im{ii) are all approx- 
imately equal, so this value can be compared directly to 
r in Eq. 2.) The two measured widths are fairly close, 
but are not equal within the experimental uncertainty. A 
possible cause of this difference will be discussed later in 
this paper. 

This Gaussian fit is strong evidence for inhomogeneous 
broadening of the linewidth. However, the cause of this 
inhomogeneous broadening is unclear. The most obvious 
source is local magnetic fields, which could be due to the 
dipolar field of neighboring clusters or the nuclear mo- 
ments of the manganese atoms. The local dipolar field 
has not, to our knowledge, been calculated because it 
depends on the random orientations along the ±z direc- 
tions of the magnetic moments of all the other clusters. 
The calculated value of the hyperfine field varies depend- 
ing on the type of coupling that is assumed between the 
electron spins, but-its maximum possible value ranges 
from 270 to 533 G,Efl and its FWHM has been estimated 
at 280-380 G.^ 

The effect of a local field would be to raise or lower 
the transition energy between the ±10 and ±9 levels, 
Ai^io.g. From the Hamiltonian, we see that Ai?io,9(i?) 
= A£'io,9(-B = 0) ± ^BgBz- For this material, g is very 
close to 2, so to explain a linewidth of 5.5 GHz requires 
a magnetic field distribution with width of 0.20 T. 

This local field distribution seems to be ruled out by 
the measurements of the width of the tunneling peak 
performed by Friedman et alx^ Starting with a sample 
cooled in zero applied field, a small magnetic field was 
applied parallel to the z-axis and the relaxation rate of 
the magnetization toward its equilibrium value was mea- 
sured. They measured the full width at half maximum to 
be 236 Oe at 2.6 K. If the local field had a full width at 
half maximum of 0.20 T, as implied above, then this nar- 
row peak could never have been observed. In addition, 
the line shape measured by Friedman et al. was clearly 
Lorentzian. This implies that the inhomogeneous broad- 
ening due to local magnetic fields must be significantly 
smaller than the homogenous effect of lifetime broaden- 
ing. We note that Friedman et aL's interpretation of this 
homogeneous broadening as being due to a lifetime of 250 
ps does not imply that Lorentzian tails should be visible 
in the transitions between m = ±10 and m = ±9, since 
the lifetime of states with larger |m| could be expected to 
be longer than that of the states involved in tunneling. 
However, the absence of broadening of Friedman et al.^s 
data due to local magnetic fields is puzzling. 

One possible explanation is that the width measured 
in tunneling experiments reflects collective effects. The 
relaxation rate measured by Friedman et al. is a fit to 
the exponential relaxation that is observed after an initial 
non-exponential relaxation. As was pointed out in Ref. 



[TtI , relaxation must actually be a collective process in 
which the relaxation of one molecule changes the dipole 
fields of its neighbors, thereby bringing new molecules 
into resonance. Therefore, it is possible that the width in 
the relaxation rate would be different from the instanta- 
neous width measured in photon absorption experiments. 

In fact, our calculations of the dipolar field show that 
collective effects must decrease the width of this tunnel- 
ing peak. Using the atomic positions tabulated in Ref. 
[ll] , we calculated the effective field on one cluster due to 
each of its neighbors. Because the clusters have sizes that 
are not negligible compared to their separations, this was 
done by calculating the dipolar field due to the electronic 
spin of every Mn atom in the neighboring cluster with 
moment nij = gfiBSjZ, where we choose Sj = ±2 for the 
eight Mn'^^ ions and Sj = —3/2 for the four Mn^+ ions to 
refiect their antiferromagnetic coupling. The fields due 
to all twelve Mn ions located at positions rj in the neigh- 
boring cluster are summed at the position of every Mn 
atom in the cluster under consideration: 



B, 



^ 47r ^ rf, 
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zi-l (3) 



where is the relative position Vj — r j . (We assume that 
the electronic spins are localized on the atomic position, 
ignoring the actual electronic density.) An effective lon- 
gitudinal field is then calculated: 



eff 



"Itotal 



(4) 



where mtotai = lOg/is- 

In order to correctly calculate the standard deviation 
of the dipolar field due to N neighbors exactly, it would 
be necessary to consider all 2^ possible arrangements of 
their spins. However, we can obtain a lower limit on this 
standard deviation simply by considering the neighbor- 
ing clusters that contribute the largest fields. As shown 
in Table 1, there are only 28 neighboring clusters that 
contribute an effective field greater than 10 G. Assuming 
there is no ordering of the dipoles, these fields add ran- 
domly to give a FWHM of 520 G, as shown in Figure 3. 
The combination of this field and the hyperfine field is 
not sufficient to explain the width observed in our spec- 
troscopy. However, we note that it is considerably larger 
than the width of 236 Oe FWHM observed by Friedman 
et at, implying that collective effects narrow the observed 
tunneling width. 

One way we could understand the spectroscopy and 
the tunneling results simultaneously is if there were some 
source of inhomogeneous broadening other than local 
magnetic fields. For example, if a, the anisotropy con- 
stant in the Hamiltonian, varied slightly (~ 1%) between 
clusters due to small configurational changes of the lig- 
ands surrounding them, then this variation would be seen 
in the linewidth of the photon-induced transition, but 
not in the zero field tunneling. The photon absorption 
experiments would be sensitive to variations in a, but the 
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TABLE I: Fields due to neighboring clusters are calculated 
using Eqs. ^ and (^. All clusters that result in effective 
fields greater than 10 G are listed. The relative location is 
given in terms of the unit cell dimensions. 



location 



field 



(0,0,±1) 

(±1,0,0) 

(0,±1,0) 
(±0.5,±0.5,±0.5) 

(0,0,±2) 
(±0.5,±0.5,±1.5) 
(±1,±1,0) 



95.9 
56.9 
56.9 
28.4 
19.2 
18.6 
15.6 
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FIG. 3: Histogram of the distribution of fields resulting from 
all possible orientations of the neighboring clusters listed in 
Table 1. The line is a Gaussian with FWHM = 520 G. 



width seen in the zero field tunnehng experiments would 
be due to homogeneous broadening mechanisms, either 
due to the tunneling time itself or the interactions with 
phonons, as proposed in Ref. Calculations by Ped- 
erson et al. suggest that the arrangement of the ligands 
is important in determining the anisotropy energy.E^ We 
note that isomerism of the ligands has been shown to ex- 
ist in these molecules and results in large variations in the 



parameters in the HamiltonianjliSEI the variations in 
the ligand positions responsible for the above effect would 
have to be less dramatic. It is possible that proximity to 
the isomers could cause neighboring clusters to have some 
small variation in their Hamiltonians. We also note that 
recent work of Chudnovsky and Garanin postulates dis- 
locations as the source of spin tunpfiling and calculates 
the effect of these dislocations on a £3 The distribution in 
anisotropy constants calculated in that paper is in very 
good agreement with the distribution we observe in spec- 
troscopic measurements. It is slightly narrower than the 
distribution observed by Mukhin et al, but since the den- 
sity of defects was chosen arbitrarily, this could simply 
imply that the defect density was higher in the sample 
measured by Mukhin et al. 

Such a distribution of a would have no effect on the 
width of the relaxation peak in zero field, since in zero 
field all levels ±m are degenerate regardless of a. How- 
ever, tunneling peaks at non-zero fields would be broad- 
ened, since (neglecting 0) the field iJ„ = —an/gjiB that 
brings the levels m and —m + n into resonance is propor- 
tional to a. It would be difficult to observe this broaden- 
ing in tunneling measurements, since the term (38^ brings 
different levels into resonance at different fields, so that 
if more than one level m is involved in the thermally- 
assisted tunneling, the relaxation peak will be broadened. 
However, we note that the width of the relaxation peaks 
in non-zero fields observed in Ref. is sufhciently large 
to accommodate the required distribution of a. 

In conclusion, we measure the linewidth of the sin- 
gle photon transition from m = ±10 to to = ±9. This 
linewidth, while surprisingly large, is in agreement with 
that observed by Mukhin et ala The combinationjjaif this 
linewidth and that measured by Friedman et alx3 sug- 
gests that the anisotropy constant varies between clus- 
ters. 
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